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The dynamics of entanglement and quantum discord between two identical qubits strongly interacting with a 
common single mode leaky cavity field have been investigated beyond the rotating wave approximation (RWA) 
by using recently derived Lindblad type quantum optical master equation [F. Beaudoin, J.M. Gambetta, A. Blais, 
Phys. Rev. A 84, 043832 (201 1)] that can describe the losses of the cavity field in a strong atom-field coupling 
regime. Contrary to previous investigations of the same model in the dissipative regime by using the standard 
Lindblad quantum optical master equation in a strong-coupling regime, the atom-field steady states are found 
to be cavity decay rate independent and have a very simple structure determined solely by the overlap of initial 
atomic state with the subradiant state which is valid for all coupling regimes. Non-RWA dynamics are found 
to have remarkable effects on the steady state quantum discord and entanglement that cannot be achieved under 
RWA conditions, for instance, they can induce steady state entanglement even for the initial states that have 
no overlap with the subradiant state. Moreover, the non-RWA dynamics are found to reverse the initial state 
dependence of steady state entanglement and quantum discord contrary to the RWA case. 
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I. INTRODUCTION 

The description of light-matter interaction is the fundamen- 
tal question investigated in circuit and cavity quantum elec- 
trodynamics (QED). The single-mode spin-boson model, the 
so-called quantum Rabi model, is the simplest possible physi- 
cal model that describes the interaction of a two-level atom (a 
qubit) with a quantized single mode electromagnetic field (a 
harmonic oscillator) [Q]]. Indeed, this simple model has been 
used to understand a wide variety of phenomena in quantum 
optics and condensed matter systems, such as quantum dots, 
trapped ions, superconducting (artificial) qubits, optical and 
microwave cavity QED, among others 12|-|4|]. Despite its old 
age and numerous investigations over the last few decades, the 
spectrum of the Rabi model has been given only recently in a 
closed form as the solution of a transcendental equation with 
a single variable Q5|] . On the other hand, the Rabi model is 
quite simplified and can be solved analytically when the ex- 
citation non-conserving terms, the so-called counter-rotating 
terms that are responsible for simultaneously exciting or de- 
exciting the atom and field, are ignored Jg]. This approxi- 
mation is known as the rotating wave approximation (RWA) 
which reduces the Rabi model to the Jaynes-Cummings model 
and has been widely used in quantum optical settings. The 
RWA is valid only at weak atom-field coupling strengths and 
nearly resonant conditions. Indeed, the typical quantum optics 
experiments are performed at microwave to visible light range 
of electromagnetic wave spectrum where the atom-field cou- 
pling strength is several orders of magnitude smaller than the 
transition frequencies of the atom and the field, so that RWA 
is inherently justified for those parameters. Despite being an 
approximation, the analytical achievements of RWA have suc- 
cessfully explained many novel phenomena and experimen- 
tal results, such as Rabi oscillations 0], squeezing Ht], non- 
classical |0] and Fock fioll states, collapse and revivals of 
atomic inversion ifTTll . and entanglement between atoms as 
well as atom and field systems 111211 . However, with the re- 
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cent developments in the area of circuit and cavity QED sys- 
tems 0,01, ultra and deep strong light-matter couplings be- 
came experimentally achievable which makes it necessary to 
use full the Rabi model to explain all the dynamical and sta- 
tistical properties of these systems. 

The discussion of the role of counter-rotating terms on 
the atomic dy namics has received plenty of attention, re- 
cently itO - tlal . Although the effects of these terms are ex- 
tremely small and averages to zero over a very short time scale 
for the weak coupling regime, they, nonetheless, cause signif- 
icant shifts in the dynamics of the system at intermediate to 
deep interaction strength region. Generation of photons rfl3ll 
and entanglement lfl4tl from zero excitation initial state, bifur- 
cation in the phase space lfl5tl . a fine structure in the optical 
Stern-Gerlach effec t iflfjll . Bloch-Siegert shift Jnl and chaos 
in the Dicke model flal are the remarkable examples of these 
novel effects and none of them are observed under RWA. On 
the other hand, the significance of the counter-rotating terms, 
especially for multipartite and dissipative systems, is still un- 
der debate and needs further investigation. 

Any realistic quantum system inevitably interacts with the 
outside world which leads to decoherence. In circuit and cav- 
ity QED systems with atom-field interactions, this unavoid- 
able interaction leads to qubit relaxation and dephasing as well 
as cavity relaxation. The standard Lindblad quantum optical 
master equation JUt] is the general form of a Markovian mas- 
ter equation that can describe the irreversible and non-unitary 
processes of the atom-field system {see Eq. (14) in Ref. l20ll }. 
While deriving the standard dissipators in this master equa- 
tion, the qubit-field interaction is neglected and therefore sep- 
arate interactions of the involved parts with the environment 
are taken into account. For the weak atom-field coupling 
regime where RWA can be safely applied, this master equation 
works well and can be used to predict the experimental results 
in cavity and circuit QED systems J2lll . On the other hand, it 
is controversial to use the standard Lindblad master equation 
without any formal proof to describe the dephasing and relax- 
ation processes in strong coupling regime where RWA is no 
longer valid. In fact, qubit-field coupling becomes influential 
and leads to unphysical effects 12011 : even at zero tempera- 
ture, the standard master equation with the Rabi model pro- 
duces spurious qubit flipping and creates excess photons in 
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the atom-field system. It was shown in Ref. Ir20ll that it is pos- 
sible to correct the predictions of the Rabi model in the dissi- 
pative case by carefully deriving a new Lindblad type dissipa- 
tor which involves transitions between the eigenstates of the 
Rabi model. The improved Lindblad master equation given 
in Ref. lEoTl is expected to be a new avenue in the studies of 
strong light-matter interactions including dissipative and mul- 
tipartite systems and to help understanding the role of counter- 
rotating terms under dissipation and dephasing processes. 

In the present study, we extend the above ideas by study- 
ing a model of two cold-trapped atoms resonantly interacting 
with a single-mode leaky cavity field. This model is quite 
versatile and can be directly verifiable in circuit and cavity 
QED setups. Under RWA, the dynamics of the atom-field 
system under losses have been investigated thoroughly and 
the atom-field steady states have a very interesting structure 
determined solely by the overla p be tween the initial state and 
the decoherence free state II22h25I1 . On the other hand, the 
same model has also been investigated under non-RWA con- 
ditions by using the standard Lindblad master equation ll26l - 
l29ll . Since the master equation with standard dissipators can- 
not explain properly the effect of dissipation in the strong cou- 
pling regime, the reported role of the counter-rotating terms 
should be reconsidered. By using the improved Lindblad mas- 
ter equation given in Ref. |f20j| and the structure of atom-field 
steady states under RWA, we provide a clear picture for the 
role of counter-rotating terms on atom-field steady states and 
also give their analytical expressions for the general class of 
two-qubit initial states, the so-called "extended" Werner-like 
states 1 30 1 which may reduce to the Werner-like mixed states 
or to the Bell-like pure states for certain conditions. We also 
reexamined some of the main results obtained by the standard 
Lindblad master equation in the strong atom-field coupling 
regime and revised them by using the results obtained by the 
improved Lindblad master equation. Beside this, we also an- 
alyzed the role of counter-rotating terms on the dynamics of 
quantum correlations, such as entanglement and quantum dis- 
cord between two atoms and found that the counter-rotating 
terms can, peculiarly, reverse the initial state dependence of 
the steady state entanglement and quantum discord contrary 
to the RWA case. For RWA dynamics, the overlap of ini- 
tial state with the decoherence free state is crucial to obtain 
a long-lived quantum correlated state under losses 1I22U25I1 . 
We showed that the counter-rotating term contributions break 
down this description and lead to long-lived quantum corre- 
lated state without this overlap. 

The paper is structured as follow. In Sec. |II] we present the 
system Hamiltonians for both RWA and non-RWA conditions 
and the master equations that describe the losses of the cavity 
field. The correlation measures, quantum discord and entan- 
glement, are also briefly discussed in this section. In Sec.lHTl 
we present our results. We conclude with a brief summary of 
important results in Sec. |IV] 



II. HAMILTONIAN AND THE MASTER EQUATION 

The Rabi model describes the interaction of a two-level 
atom with a single mode electromagnetic field. The Hamil- 



tonian can be given as (with H = 1 throughout) |Q]]: 

H R = ~^ a z + UrCi^a + g(a + + cr_)(a + at), (1) 

where ui a (u r ) is the transition frequency for the atom (cavity 
field), a z is the qubit energy operator, a± are the spin-flip 
operators and g is the atom-field coupling constant. 

The terms aV+ and acr_ that appear in the interaction part 
of the Rabi Hamiltonian are called counter-rotating terms and 
they describe "virtual processes"; aer_ describes a process in 
which a photon is annihilated in the cavity mode as the atom 
makes a transition from its excited to the ground state, while 
a'<T+ describes the creation of a photon in the cavity mode as 
the atom makes an upward transition. So, the counter-rotating 
terms do not conserve the total number of excitations in the 
system. 

The virtual processes have observable effects when g is 
comparable with us a and lu t , i.e., in strong atom-field cou- 
pling regime, however when the weak coupling, g << uj a , ui r , 
and the nearly resonant, |A| — \u> a — cu r \ << uj a , cj r , condi- 
tions are simultaneously satisfied, these effects become negli- 
gible and the counter-rotating terms can be neglected in the 
interaction part of the Rabi Hamiltonian. This approxima- 
tion is known as the rotating wave approximation (RWA) and 
the Rabi model reduces to the well known Jaynes-Cummings 
model given by the Hamiltonian @]: 

Hjc = -^-<J Z + u r a) a + g{o + a + (j-o)). (2) 

In opposition to the Rabi Hamiltonian, the total number of 
excitations in Hjc is a conserved quantity which makes Hjc 
analytically accessible. 

In cavity QED systems, there would be photon loss from 
the cavity due to the imperfections in the cavity mirrors. At 
absolute zero temperature (T = 0) and under Markov ap- 
proximation, this can be represented by the standard Lindblad 
quantum optical master equation lfl9ll : 

p = -i[H,p] + KD[a]p, (3) 

where D[m\p = 5 (2mpm^ — pm^m — m^mpj, k is the 
photon leakage rate, m is the appropriate dissipator opera- 
tor, and p is the atom-field density matrix. The first term in 
Eq. ® describes the unitary Hamiltonian dynamics, while the 
second part accounts for the dissipation. At the weak-coupling 
regime, where RWA can be applied, the dissipator part of the 
master equation is well known and involves field annihilation 
operator, a, which is responsible for the leakage of photons 
from the cavity. Indeed, at weak-coupling regime where RWA 
can be applied, this master equation can be safely used and the 
dissipation can bring the system to the ground state, \g0), of 
Hjc- 

Although it has been widely investigated by many groups, 
recently J3I11MIIMI, 

using the above dissipator in the 
dissipative Rabi model in an ad hoc manner has been shown 
to lead to unphysical consequences [20]; at T = 0, the dissi- 
pation will generate excess excitations when even no energy 
is added to the system and consequently the dissipation will 
drive the system out of the ground state of Hr. In fact, the 
main failure of this master equation in the strong-coupling 
regime is to neglect the qubit-field coupling when deriving 
the dissipator part of the master equation. In Ref. ll2o"ll . the 
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qubit- field coupling have been included in the derivation of 
the master equation and the new Lindblad master equation that 
describes the losses of the cavity field at T — and strong 
coupling regime under Markov approximation can be given 
as: 



P = -i[H R ,p] 



53 T{ k D[\j) (k\]p, 

J,k>j 



(4) 



where are the relaxation coefficients and are related to 
the spectral de nsity of the bath and the system-bath cou- 
pling constant 02011 . It can be simplified as (H, T£ fc = 



U\(a 



\k) | 2 . Here \j) and ujj are the eigen- 



states (dressed basis) and the eigenvalues of the Rabi Hamil- 



tonian, respectively, i.e., Hr \ j) 



\j) and the eigenstates 



are labeled according to increasing energy, i.e., label \j) such 
that ujk > ojj for k > j. Indeed, as shown in Ref. 12011 . the 
new dissipator drives the atom-field system to the true ground 
state. Also, note that the standard dissipator of Eq. (f3} can be 
obtained from that of Eq. © in the limit g — > 0. 

Since the strong atom-field coupling regime became rele- 
vant with the development of the recent technology (HI 01 and 
a number of investigations of the dissipative Rabi model with 
standard dissipator have been re port ed in the literature in the 
strong-coupling regime ll2^l - l2^l3Tij34ll . a reexamination of 
these studies in the light of the improved master equation <j4j 
is due. Toward that goal, in the present study, we will analyze 
the dynamics of entanglement and quantum discord between 
two identical qubits (A and B) interacting with a common sin- 
gle mode leaky cavity field. The Hamiltonians for non-RWA 
and RWA cases can be found by just replacing a z — > aj- + 
and a± — > a± + a± in Eqs. (Q~|i and (fJJ, respectively. We will 
solve the master equations (fJJ and (0]i numerically for the con- 
sidered Hamiltonians and the initial states where the atoms are 
one of the four extended Werner-like states ll30ll and the cavity 
field is in its vacuum: 



P*±(0) 
/>*±(0) 



1-r 
4 

1 - r 



lAB+r\*t)(*t\ 



®|0) (0|, 

® |0> (0| , (5) 



where Tab is the 4x4 identity, r (0 < r < 1) is the purity of 



the initial states, I $ 



±\ 



\eg)± Vl - a 2 \ge) and \^>±) 



a\ee) ± \/l — a 2 \gg) are the Bell-like states, and a (0 < 
a < 1) is called the degree of correlations. The extended 
Werner-like states are mixed and reduce to the well known 
Werner-like mixed states at a = l/y/2 and to the Bell-like 
pure states at r = 1. The atomic reduced density matrix can be 
obtained by taking a partial trace of atom-field density matrix 
over the cavity degrees of freedom. The above initial states in 
the standard basis {|1) = \ee) , |2) = \eg) , |3) = \ge) , |4) = 
\gg)} have an X structure in the atomic space with non-zero 
elements only in its main- and anti-diagonals: 



PAB = 



>11 pi4 

P22 P23 

P32 P33 

. P41 p 4 4 



(6) 



It was shown that the standard Lindblad master equation (fj) 
for both RWA and non-RWA Hamiltonians preserves the X 



structure of the density matrix ll26[ l28l l29ll . We have checked 
that the X structure of the atomic reduced density matrix re- 
mains intact also under the improved master equation (|4j. 

Since the excitation number in the system is not conserved 
for the Rabi model, it is worth to explain briefly the numerical 
technique to solve the master equations (0 and (|4]i for the two- 
qubit Rabi Hamiltonian. In the present study, we have care- 
fully checked the convergence of computed properties versus 
the dimension of the cavity Fock field. We have considered 
basis vectors of type \i,j,n), where i,j = e, g and n = 
0, 1, 2, ... M to obtain a set of differential equations for the 
density matrix elements. It is found that M rj 50 gives well- 
converged results for the highest interaction strength consid- 
ered in the present work. Moreover, depending on the inter- 
action strength, 100-150 dressed basis of the two-qubit Rabi 
Hamiltonian are used to calculate the time evolution. Dur- 
ing the numerical calculations, the basic properties of density 
matrix, such as positivity, hermiticity, and trace preservation 
properties are all monitored. 

Our next step is to determine the correlations between the 
atoms. The entanglement, a kind of quantum correlation, 
determines whether or not a given bipartite state is separa- 
ble. Entangled states are broadly accepted as a necessary 
resource for a set of quantum tasks in quantum information 
theory I' 36. [37ll . such as quantum key distribution and tele- 
portation. Entanglement can be calculated through entangle- 
ment of formation (EoF) which is the function of the entangle- 
ment monotone, the so-called concurrence, for bipartite states 
as (H: 

Ceof(pab) = -rjlog 2 ri- (1 -??)log 2 (l -?/), (7) 



where r) = 1/2 (l + Vl - C 2 ) with C = 2max{0, \p u \ - 
\JP22P33-, I/O23 1 — y/PuPii} being the concurrence for X states. 
On the other hand, quantum discord (QD) emerged as a new 
fundamental type of quantum correlation beyond entangle- 
ment (39|-|4l|]. QD as a fundamental resource is shown to 
be useful in practical quantum information tasks, such as 
DQC1 (deterministic quantum computation with one quan- 
tum bit) ll42ll. G rover search ll43ll algorithms and remote state 
preparation [44]. Its definition is based on the difference be- 
tween quantum versions of two classically equivalent defini- 
tions of mutual information. Non-zero QD signifies that it is 
impossible to extract all information about one subsystem by 
performing a set of measurements on the other subsystem. QD 
can also be calculated, analytically, for X states as ll45ll : 



Cqd(pab) = min{<5i,Q 2 }, 



(8) 



where Q 3 = h[p n + p 33 ] + J2 k =i X k log 2 A fc + Dj with A fc 
being the eigenvalues of pab and h[x] = —x\og 2 x — (1 — 
x) log 2 (l — x) is the binary entropy. Here D\ = h[r], where 

t = (l + Vll - 2(p 33 + PAi)] 2 + 4(M + |p 23 |) 2 ) /2 and 

D 2 = - J2t=i Pkk log 2 Pkk - h[p n + p 33 ]. QD and EoF are 
equal for pure states, but the relation is complicated for the 
mixed states; there are separable mixed states with non-zero 
QD l46ll . Recently, the comparison attempts for the dynamics 
of quantum discord and entanglement have received a great 
deal of attention for open quantum systems I47l - l57ll . The 
studies show that QD is much more robust compared to en- 
tanglement under dissipative and dephasing processes where 
entanglement can suffer sudden death. 
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In the following, our main concerns will be threefold. Since 
numerous studies have been recently based on the standard 
master equation (01 in the strong-coupling regime 126 291 l3ll - 
l34ll . our first thought is to reexamine some of the main results 
obtained in these studies by using the new Lindblad master 
equation (|4]i. Since the standard dissipator cannot describe 
the losses in the strong atom-field coupling regime, the role 
of counter-rotating terms on the dissipative dynamics of the 
atom-field system are far from being understood. As the sec- 
ond issue, we will provide a comprehensive picture for the 
atom-field steady states and give their analytical expressions 
for the extended Werner-like initial states. Our final focus is 
to analyze the role of counter-rotating terms on the dynamics 
of entanglement and quantum discord. In the following, we 
will restrict ourselves to the resonant case, uj a = oj r = uj. 



III. RESULTS 

We start our qualitative analysis by examining the main dis- 
agreements between the two master equations in the strong- 
coupling regime. In Fig. 1, we display the effect of cavity 
decay rate, k, on the dynamics of QD and EoF for the vac- 
uum (\ggO)) initial state at g = 0.35w. The results in the 
top plots are obtained by the use of the standard master equa- 
tion [Eq. ©], while the bottom plots display the results for the 
new dissipator [Eq. (01] in the case of two-qubit Rabi Hamilto- 
nian. Since the vacuum state contains no initial excitations to 
start with and the Jaynes-Cummings Hamiltonian conserves 
the total excitation number in the atom-field system, no type 
of correlations will be induced in the time evolution under 
RWA. On the other hand, the counter-rotating terms in the 
Rabi Hamiltonian will generate virtual excitations and conse- 
quently non-zero correlations between the atoms which can be 
trapped at a non-zero value under cavity losses for both cases. 
If we study the effect of cavity losses on the dynamics of cor- 
relations by using the standard dissipator [Figs. 1(a) and 1(b)] 
as was done in Refs. l29l[3Tl - t33ll . the results seem to be quite 
surprising; as n increases, although the disentanglement time 
and the maximum of EoF decreases, the steady state QD, pe- 
culiarly, increases with n and is saturated in the maximum 
possible value (CQD max = 1/3) in the absence of entangle- 
ment for a very bad quality cavity case (k/oj > 20) 12911 . 
Although these results seem to indicate a surprising possibil- 
ity of dissipation enhancement of quantum discord rf29ll (as 
well as atomic energy as found in Refs. l3ll]33ll ) to its max- 
imally possible value for a separable state, its source is the 
error in the dissipator which creates excess excitations in the 
steady state. Considering the role of k on the dynamics of 
QD [Fig. 1(c)] and EoF [Fig. 1(d)] for the new dissipator 
demonstrates that the steady states are k independent; QD 
and EoF reaches 0.069 and 0.0086 in the long time limit, re- 
spectively, regardless of the magnitude of k. In fact, n only 
determines how fast the correlations reach their steady state 
values. This is quite understandable, since k is just a pa- 
rameter that determines how fast the photons can escape from 
the cavity mirrors, it should control the speed of reaching the 
steady state rather than its magnitude in the long time limit. 
Bath decay dependent steady state shifts of several proper- 
ties of the system in the strong-coupling regime have been 
recently investigated by using the standard dissipator ll29ll3ll - 



Moreover, Ma et al. in Ref. JH] studied the dynam- 
ics of entanglement between two qubits strongly interacting 
with a common Lorentz-broadened cavity mode at zero tem- 
perature by using a more general master equation approach, 
the so-called hierarchical equation method without employ- 
ing rotating wave, Born, and Markovian app roximations . One 
should note that the approach in Ref. Ipal is quite different 
compared to the standard Lindblad master equation case, but 
still the asymptotic values of entanglement are found to be dis- 
sipation strength dependent, which is in contradiction to our 
findings. The simple example shown in Fig. 1 indicates that 
some of the previously obtained results with the Rabi model 
need to be reexamined. As shown in Fig. 1(b) for the standard 
dissipator case, the atomic steady states become separable as 
K > 0.05w rt29tl . while the steady states can posses entangle- 
ment for the new dissipator case [Fig. 1(d)]. 

Now, we will consider the steady states of RWA and non- 
RWA dynamics. The steady state density matrix under cav- 
ity losses and RWA dynamics investigated by the standard 
master equation (0 are known to be K-independent I22I l23ll . 
As shown in the bottom plots of Fig. 1, they are also k- 
independent in the case of non-RWA dynamics with the im- 
proved dissipator. Therefore, in the following, we will set 
K = O.Ilo. 

The dynamics under RWA with losses given by Eq. (0 have 
been investigated by many groups, which indicate that the 
atom-field steady states have a simple structure. We follow 
the ideas developed in Refs. miHt] to determine those states. 
The cavity decay tends to throw out any initial excitations 
from the cavity mirror and consequently it drives the atom- 
field system to the ground state ( | <?.g0) ) of two-qubit -Hjc- 
On the other hand, there exist a highly entangled subradiant 
state (decoherence-free state) for this model due to the inter- 
action of the atoms with a common environment which leads 
to an indirect qubit-qubit interaction 1I22I, l23ll . This subradiant 
state for this model can be simply determined by finding the 
eigenvector of two-qubit Hjc for zero eigenvalue, which is 
the maximally entangled (Bell) state in the atomic subsystem 
given as |$~0) = ^= (|eg0) — |<?e0)). If an initial atomic 
state has a non-vanishing overlap with the subradiant state, 
b = ($~ I pab(0) \&~), the amount of overlap determined by 
b would be trapped in the subradiant state, while the remaining 
part, 1 — b, would decay to the ground state. Therefore, de- 
pending on the specific initial state 15^1 . the atom-field steady 
states would be a combination of \gg0) and |3>~0) which can 
be written as: 

P ss = (1 - b) \gg0) (gg0\ + b |*"0) ($"0| . (9) 

This, indeed, explains why the steady states under RWA with 
losses are k and g independent. By using Eq. (0 or by study- 
ing the solution of the master equation (0, for example by the 
method of solving first order differential equations via eigen- 
values and eigenvectors or the pseudomode approach ll23ll . we 
can determine the analytic form of the steady states for the ex- 
tended Werner-like initial states (0. For p^± (0) type initial 
states, the steady states can be found as: 

(10) 
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FIG. 1. (Color online) The effect of cavity decay rate n (a) and (c) on quantum discord and (b) and (d) entanglement versus ut for \ggQ) 
initial state, k = 0.05a; (black, solid), k = 0.5a; (red, dashed), k = 2ui (blue, dotted), and k = 20a; (green, dot-dashed) in strong coupling 
regime with g — 0.35oj. The top figures are obtained by using the standard dissipator Eq. ((3}, while the bottom plots are for the improved 
dissipator Eq. {4j. 



For r = 1, i.e., for initially pure Bell-like states, I^J), 
Pq,± (0) has no overlap with the subradiant state and all the 
excitations are lost in time evolution regardless of what a is, 
while for < r < 1, there would be a non-vanishing overlap 
term with the subradiant state determined by r and the atom- 
field system will decay to a trapping state which is superposi- 
tion of \gg0) and |$~0). On the other hand, for p $ ± (0) type 
initial states, the steady states can be determined as: 

pS 4 = (~r~ ± ra ^ 1 - a2 ) 

+ T r ay/l - a 2 ^j |$~0) (*~0| . (11) 

One should note that p $ ±(0) have always a non-vanishing 
term [except for p^+ (0) at a = -^=, r = 1] with the subradi- 
ant state, so the steady states are highly a and r dependent. 

is also the subradiant state for two-qubit Hr, which 
can be proven by noting that the Hilbert space Hq of the 
Rabi model can be separated into two uncoupled subspaces 
Hi and H 2 with H x = -±= (\eg) - \ge)) g> H cavity where the 
separation is also valid under decoherence mechanism and at 
the long time limit the state in Hi becomes the dark state, 
(\eg) — \ge)) ® |0), independent of the coupling strength. 
Since Hjc is the limiting case of Hr and the master equa- 
tion <j4j can describe the cavity losses in the strong-coupling 
regime, we conjecture that such a simple structure for the 
atom-field steady states as in Eqs. ( fTOb and (fTTb should ex- 
ist in the strong-coupling regime also. Below we will provide 



numerical evidence for our conjecture. 

We start to analyze the dynamics of correlations in the 
strong-coupling regime by using the improved master equa- 
tion for the type of initial states which have no overlap 
with the subradiant state. Thus, in Fig. 2 we display QD 
and EoF versus tut for the initial state p^+(0) with r = 1, 
g = 0.2a;, and n = O.Ilj and for several a values. In the case 
of RWA, the steady state is the vacuum which has no types 
of correlations. On the other hand, as shown in Fig. 2, the 
cavity decay in the strong-coupling regime drives the atom- 
field to a state which can carry appreciably high quantum dis- 
cord (Cq D = 0.021) and entanglement (C Eo f = 0.0014) in- 
dependent of what a is. In fact, we have observed that for 
the p^,± (0) type of initial states with r — 1, the atom-field 
steady state is the ground state of two-qubit Hr regardless 
of g. One can observe that the most important difference be- 
tween the RWA and non-RWA dynamics at the long-time limit 
is the existence of non-zero QD and EoF created by the virtual 
processes; for RWA dynamics the necessary condition for the 
long-lived correlated state is the overlap of initial state with 
the subradiant state 1I22I - I25I1 . while this can be achieved with- 
out the overlap under non-RWA conditions because of the en- 
tangled nature of the ground state. In the insets of Fig. 2 we 
display the g dependence of the ground state QD and EoF. 
Both QD and EoF as a function of g/u display a peculiar 
structure; they have a single maximum at the intermediate 
values of g/u. Furthermore, g/u> dependence of EoF is al- 
most Gaussian. For small g/oj, where the RWA regime is 
valid, the QD and EoF are nearly zero, since the ground state 
is close to the vacuum state, while the correlations are high 
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FIG. 2. (Color online) Effects of degree of correlations a on (a) QD 
and (b) EoF versus urt for p^+(0) initial state with r = 1, g = 
0.2a;, At = 0.1a;, and a = (black, solid), a = 0.2 (red, dashed), 
a = l/y/2 (blue, dotted), a = 0.8 (green, dot-dashed) and a = 
1 (yellow, dot-dot-dashed). Here the results are obtained by using 
two-qubit Hr and the master equation, Eq. ((4). Inset: QD and EoF 
of the ground state versus g/u). Around g/u fts 1.3, the ground 
state of the Rabi model becomes almost twofold degenerate. So, we 
consider only g/u £ [0, 1.2] range in the insets. 



FIG. 3. (Color online) Steady state (a) QD and (b) EoF versus 
initial purity r for p^±(0) initial states under RWA (red, dashed) 
and non-RWA dynamics with g = 0.05a; (yellow, dot-dot-dashed), 
g — 0.2a; (black, solid), g = 0.35a; (blue, dotted), g = 0.6a; (green, 
dot-dashed) and g — 0.75a; (black, long-dotted). Note that for 
g < 0.01a;, the results obtained by the new master equation ((4} co- 
incide with the RWA case (red, dashed) given by Eq. d 1 0b and the 
steady states for both cases are a independent. 



between 0.2 < g/u < 0.7. Note that for the g/u regions, 
where EoF is absent, non-zero QD still exists especially for 
g/oj > 0.8. Also note that the ground state QD and EoF 
show peaks at different g/u values; QD is maximum nearly at 
g/u> ps 0.53, while EoF is at g/u ps 0.48. On the other hand, 
the decrease in both QD and EoF as a function of g /u at very 
strong-coupling regime might be an indication of the inability 
of the Rabi model to account for the possible nonlinear effects 
at large g values, which may require higher order terms in the 
description of atom-field interactions Il60ll . 

Special attention should also be given to the initial Bell 

= -75 (N) + |<7e)) and Bell-like |tf±) = a\ee) ± 



VI — a 2 \gg) states in the strong-coupling regime, since the 
entanglement dynamics have been extensively analyzed by 
using the standard dissipator (0 under losses N26h29I1 . For 
k > O.lu, entanglement between atoms becomes zero or use- 
less (smaller than 10~ 4 ) for a wide range of g, especially be- 
tween 0.2 < g/u < 1 under the standard dissipator case ll26l - 
l29tl . Contrary to those findings, one can see that such states 
can contain appreciably high EoF for 0.2 < g/u < 0.7 inde- 
pendent of k as shown in the inset of Fig. 2(b), when one uses 
the improved master equation. 



Now we consider the dynamics of initial states that have 
non-zero overlap with the dark state, |$~0). The effect of r 
for p^± (0) initial states on the dynamics of QD and EoF un- 
der losses have been investigated and the r dependence of the 
steady state correlations are plotted in Fig. 3 for both RWA 
and non-RWA cases. Before discussing the role of counter- 
rotating terms on the steady state quantum correlations as a 
function of r, we should stress that a detailed analysis demon- 
strates that the atom-field steady states in the strong-coupling 
regime are in the form of Eq. (ITOi with just \gg0) replaced by 



the ground state, 
as: 



ggO), of two-qubit Hr and can be written 



ggO)(ggO 



|$"0) <$"0| 



(12) 



Since the dark state of two-qubit Hjc is also the dark state 
of two-qubit Hr, the form of the atom-field steady state for 
non-RWA Hamiltonian [Eq. ( fT2b l is in expected form based 
on the same arguments that lead to Eq. ([Tot . More precisely, 
the overlap of the atomic initial state with the subradiant state 
is trapped in the subradiant state, while the remaining part 
decays to the ground state of the open quantum system un- 
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der consideration. Also note that in the small g limit ggO 

reduces to \ggO). Therefore, the overall role of the counter- 
rotating term contributions in the steady states stems from 
the ground state of the two-qubit Rabi Hamiltonian. We have 
checked the form of Eq. (fT2l i extensively by numerical means 
and found that it holds for all the initial states considered in 
the present work. 

The initial purity dependence of steady state QD and EoF 
displayed in Fig. 3 for the RWA case shows an interesting fea- 
ture; steady state values of both correlation measures show 
an inverse dependence on r. Since r is a measure of initial 
mixedness of the state and the dissipation dynamics lead to 
an increase of the mixedness of the initial state, this finding 
might seem to be contradictory, but it stems from the fact that 
as r decreases the overlap between the initial state and the 
maximally correlated dark state increases and the final steady 
state would contain a larger fraction of |$~) which leads to 
higher QD and EoF. In the strong-coupling regime, the ground 
state is g dependent, which makes QD and EoF dependent 
on g also. In fact, the counter-rotating contributions compli- 
cate the simple RWA picture and can induce separable atomic 
steady states for a wide range of r. Similar to the RWA case, 
the steady state correlations can display inverse r dependence 
for a wide range, but contrary to the RWA case, the correla- 
tions can also be linearly r dependent in the strong-coupling 
regime. Moreover, a generic result can be obtained from the 
comparison of the magnitudes of QD and EoF in Fig. 3; for the 
regions where EoF is zero (especially for strong couplings), 
appreciably high QD can still exist in the steady states. 

To further elucidate the role of initial state and virtual pro- 
cesses on the steady state quantum correlations, we display 
QD versus a for (0) with r = 1 in Fig. 4 and QD ver- 
sus r for the same initial states with a = 0.2 in the insets 
of Fig. 4. Similar to the p <S/ ± (0) case, our detailed numerical 
analysis shows that the atom-field steady states under non- 
RWA dynamics for p^± (0) type initial states are in the form 
of Eq. (fTTT i with the vacuum state replaced by the ground state 
of the Rabi Hamiltonian, i.e., 



(a)p* ; (0) 2\ 

k RWA . 20 . 

V g=0.05w a 0i5 _ 

v g=0.75w 010 




0.0 



0.0 0.2 



0.0 0.2 0.4 0.6 0.8 1.0 



0.4 0.6 

a 



1.0 



FIG. 4. (Color online) Steady state QD versus degree of correlations 
a for (a) p s +(0) and (b) p s -(0) initial states with r = 1 under 
RWA (red, dashed) and non-RWA dynamics with g = 0.05oj (yel- 
low, dot-dot-dashed), g — 0.2oj (black, solid), g = 0.35oj (blue, 
dotted), g — 0.6oj (green, dot-dashed), and g = 0.75u; (black, long- 
dotted). Inset: Steady state QD versus r for (a) p s + (0) and (b) 
p s _ (0) initial states with a = 0.2 and the same g values used in 
the subfigures. Note that for g < O.OIoj, the results obtained by the 
new master equation 101 coincide with the RWA case (red, dashed) 
given by Eq. i fTTV 



IV. CONCLUSION 



n ss 



3-r 
4 

1 + r 



± ra\/ 1 — a 2 



=F ra 



ggO)(ggO 



|$-0)($-0|. (13) 



Concisely analyzing the steady state QD dependence on the 
initial state parameters under RWA (red, dashed line) signifies 
that the magnitude of QD directly depends on the amount of 
overlap with the subradiant state. On the other hand, strong 
virtual processes can stir the simple RWA picture and inter- 
estingly can reverse the initial state parameter dependence of 
the steady QD compared to the RWA case. The a and r de- 
pendence of steady state EoF is qualitatively similar to that 
of QD with a difference as previously indicated in Fig. 3(b); 
EoF can suffer death in the steady states for a wide range of a 
and r under non-RWA dynamics, especially at g = 0.6w and 
g = 0.75a;. Thus, they are not plotted here. 



We have studied a model of two cold-trapped two-level 
atoms strongly coupled to a single-mode leaky cavity field 
and initially prepared in extended Werner-like states. The dy- 
namics of the atom-field system and quantum correlations be- 
tween atoms, such as quantum discord and entanglement in 
the strong-coupling regime, have been investigated by using 
the recently derived Lindblad type master equation which of- 
fers quantumjumps among the eigenstates of the Hamiltonian 
of the open quantum system under consideration. Our results 
demonstrate that for extended Werner-like initial states, the 
amount of overlap between initial and dark states remains in- 
tact in the time evolution, while the other part decays to the 
ground state of the atom-field system. This result is found 
to hold for all atom-field coupling strengths and coincides at 
the weak-coupling regime with the results under RWA stud- 
ied by the standard Lindblad master equation. The presence of 
counter-rotating terms are found to break down some of the re- 
cently predicted results for RWA dynamics, for example, they 
lead to quantum-correlated atomic steady states without the 
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overlap of an initial state with the dark state which is crucial 
for RWA conditions. Moreover, the virtual processes can re- 
verse the initial state parameter dependence of steady state QD 
and EoF contrary to the RWA case which is strongly initial- 
dark state overlap dependent. However, we should stress here 
that the strong coupling regime seems to have quantum advan- 
tages compared to the RWA case only for initial states which 
lead to low (Cqd < 0.05) and zero discordant atomic steady 
states for the RWA case. Specifically, the strong virtual pro- 
cesses can induce separable atomic steady states. 

The results obtained in the present study are quite simple 
and general. More importantly, the developed ideas can be ex- 



panded for local or global environment cases including multi- 
partite and other decoherence channels, such as dephasing and 
spontaneous emissions of qubits and fields. For instance, for 
a model system with two qubits locally trapped inside leaky 
cavity fields l6lll . the dark state no longer exists. Therefore, 
it is expected that the atom-field steady states under cavity 
decay would always be the ground state of the system Hamil- 
tonian. However, a more detailed discussion for this subject 
goes beyond the aim of the present study. 
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